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Abstract A method is developed for semi-automated 
classification of SAB, images of the tropical forest. Infor- 
mation is extracted using the wavelet transform (WT). The 
transform allows for extraction of structural information in 
the image as a function of scale. In order to classify the 
SAR image, a Desicion TVee Classifier is used. The method 
of pruning is used to optimize classification rate versus tree 
size. The results give explicit insight on the type of infor- 
mation useful for a given class. 
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1. Introduction 

The work presented is part of the Global Rain Forest 
Mapping project (GRFM). Recent work have shown the 
potential and use of Radar remote sensing for monitoring 
of tropical forests (l] f [2], The objective of our team is 
to produce an automated classification algorithm used for 
automated mapping and monitoring of tropical forest on 
a global scale using SAR. We present a method for semi- 
automated classification of SAR images. 

The set of data covers the entire tropical belt with JERS- 
1 SAR data acquired at two dates. It was provided by 
NASDA (japanese space agency). The resolution of the 
data was reduced to 100 m before mosaicking. In this paper, 
we will focus on a part of Gabon where ground truth is 
available. 

Various methods exist to classify images using training 
areas. One of the most popular is the “Maximum Likeli- 
hood” method (ML). Other methods such as neural net- 
works and decision trees offer advantages related to the un- 
derlying assumptions which are not needed in the last two 
methods as well as the explicit information extraction on 
the various data contributions. In this paper, the decision 
tree method was used [3]. The tree uses training samples 
to grow, based on node impurity measure. A split is chosen 
f it maximizes the the reduction of impurity (*(t)) of the 
parent node t with respect to the impurity of the splits. 

Ai(0 = i{t) - Pl i(t t ) - prl(tr) (1) 

Pl and p r are the proportion of samples in node t that reach 
into child nodes l and r. 

The Gini criterion was chosen, as it was shown to be 
generally efficient [3]. It is defined as: 
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It uses a rule that assigns a sample randomly selected 
from the node t to class j with probability p(j|t). A split is 
chosen to maximize Ai(f)- The probability is: 


p(j\l) = 
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where Pj is the prior probability. 

However, the tree grows indefinitely until all classes are 
separated which can lead to a very large tree. Instead, we 
will grow a reasonable tree (7 levels), followed by “prun- 
ing” of the tree. The pruning method consists in cutting 
branches wich are inefficient. The decision is based on the 
misclassification cost, R(t) = r(t)p(t), with 

r(t) = 1 - max jP (j\t) ( 4 ) 

where p(t) is the probability that any case fall in node t. 

Following the misclassification cost, a complexity func- 
tion is defined [3] that takes into account the number of 
terminal nodes nT . The weakest non-terminal node is the 
one with minumum 


9(t) = 


R(t)-R(T t ) 
nT - 1 


( 5 ) 


The original tree is iteratively cut to produce a sequence 
of subtrees T k with different complexity. The optimum-size 
tree is chosen using an independant sample which classes 
are known. This sample is input into the tree sequence and 
estimates R(T k ) are produced. The optimum tree is the one 
with minimum estimate R(T k ). 

Multiscale information from the SAR image is extracted 
using the Wavelet Transform (WT) [4]. The Transform 
produces a multiscale representation of the original images 
as well as 3 (orientation) detail images. The detail images, 
which contain the wavelet coefficients, can be considered as 
band-pass representation whereas each of the detail image 
contains information at a finite scale interval. The wavelet 
coefficients (c) are the components of the original image at 
a given scale (a) on the wavelet base (tf): 


Co =< /,* > . (6) 

Thus, structures in the SAR image are separated with 
respect to their scale. Texture maps are built from the de- 
tail images by quadratic addition of the 3 orientation detail 
images and normalisation by the low resolution local am- 
plitude. The normalisation is computed in order to correct 
for multiplicative noise contribution [5]. The low-resolution 
approximation, which is an output of Mallat’s algorithm 
contains and estimate of the local reflectivity which is used 
to estimate noise level. Normalisation was also shown to 
enhance edges [6]. 


II. Application of the method 

Mallat's WT algorithm was applied to the image shown 
on Fig. la. The texture maps were built by adding 

quadratically the detail images and normalizing by the 
low-resolution approximation resulting from the algorithm. 
Fig. lb is the texture maps at scale 200 m. As the scale 
increases, the discontinuities in the signal are enhanced. 
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Fig. L a) SAR image of part of Gabon. Training sites are shown 
by arrows, b) Texture map of the region at scale 200 m. 


Two dates were input into the classifier as well as the 
texture maps for the low water season. Fig. 2 shows the 
linear histogram stretch applied to the data to cover the 
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Fig. 2. Histograms of the input images: Low water amplitude (Full 
line), textures (dashed lines). 


(errors: Low land forest (10%), medium topography for- 
est (23%), highland forest (61%), raphia (37%), inundated 
(17%), water area (29% mixed with savanna), river (15%), 
city (24%), savanna (50 % with water), savanna temporal 
change (50%), swamp (53%), regrowth (classified as forest), 
low mangrove (mostly classified as forest 90%), high man- 
grove (50%), open forest (5%)). Grouping of classes could 
significantly reduce errors. 


III. Conclusion 

The results explicitly demonstrate the contribution of 
each input to distinction between the classes. It also shows 
the limitations of using JERS-1 SAR configuration for mon- 
itoring of tropical forests. That is, it is limited to biomass 
extremes (ex: savanna-forest) and vegetation moisture (e.g. 
inundated areas). The decision tree classifier can be used 
with different inputs. Thus, studies with different texture 
measures will be performed. We also need to compare this 
classifier to classical ones such as the ML method. 

The result errors and artefact demonstrate the impor- 
tance of choosing the training sites. This is readily seen 
when comparing misclassification Thus, these result must 
be considered when choosing the training areas, misclassi- 
fication high land forest as edges (or open forest). This is 
due to very large texture. 
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Fig. 3. Classified map of Gabon. The end nodes resulting from the 
algorithm have been grouped (in colors) to the class with largest 
proportion. Water (class 0), low land forest (class ), medium to- 
pography forest (class ) High land forest (class ), raphia (class ), 
inundated vegetation (class ), savanna (class ), Temporal change 
savanna (class ). 
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Fig. 4. Tree grown using training sample overall probabilities and 
pruning using one-class probabilities. At each node, a value cor- 
repondaing to the best split s threshold is given with the input 
parameter i (s*). The parameters are the amplitudes for low 
water (0), high water (l)\ texture 100 m (2), 200 m (3) and 400 
m (4). 



